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^ ^ ; I. INTRODUCTION 

m _ 

, It is now more than a decade since the first polarised DIS experiments ^ discovered the strong breaking of a SU (6) 
' quark model based sum-rule ||^ , and precipitated the "proton spin crisis" . Since then many polarised deep- inelastic 
scattering (DIS) experiments have reported measurements of the virtual photon asymmetry 
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on different targets (the structure functions gi and Fi are defined later). A polarised proton collider at RHIC will 
f-^ ^ soon begin to constrain the unknown polarised parton distributions even more strongly. Current and future interest 
I in this topic stems partly from the history of the "spin crisis" . 

However, polarised parton densities are also interesting because of the role they might play in future polarised 
hadron collider searches for completions of the standard model. Essentially, a large variety of physics beyond the 
standard model plays with chirality. Some of this freedom can easily be curtailed by polarised scattering experiments, 
if the polarised parton densities arc known with precision. We expect that by the end of the polarised- RHIC program 
this goal should be reached. 
; ^ ' The longitudinally polarised structure function gi is defined by 




gi{x, Q')^Cg®-y^ ei [qf + qA + 7; \ } ^ ei \Cg^g, (2) 



which is a Mellin convolution of the quark (<7/), anti-quark (g^) and gluon (g) longitudinally polarised distributions 

with the coefficient functions Cq^g. The index / denotes flavour, and Cf is the charge carried by the quark. The 
unpolarised structure function is given by a similar formula in terms of the corresponding unpolarised densities and 
coefficient functions. 

These coefficient functions and the splitting functions (which determine the evolution of the densities) are com- 
putable order by order in perturbative QCD. The former are crucial for the Bjorken sum rule [Q, connecting the 
proton and neutron structure functions, gl and g", to the neutron /3-decay constant gA, and are known to NNLO [|. 
This makes it possible to use this sum rule for precision measurement of the QCD scale The polarised splitting 
functions are known only at NLO 

In this paper we analyse the currently available inclusive DIS data in QCD and extract polarised parton distributions 
from them. In this respect our work is similar to that of 0. However, our analysis differs in several ways. For one, 
some of the data we use is more recent than the older fits. More importantly, we relax some of the assumptions which 
needed to be made in analysing the older data. We allow for flavour asymmetry in the polarised sea quark densities. 
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and let the first moment of the gluon density vary freely in the fit. Furthermore we make a detailed investigation of 
the uncertainties in these polarised gluon and sea quark densities. 

The uncertainty in gluon densities may seem puzzling in view of the fact that the Q^-dependence of the structure 
function gi involves the gluon strongly. In fact, at LO we already have — 



dgi as 



aiogQ2 2tt 



Pqq ffl + Pqg < 



(3) 



where Pqq and Pqg are polarised splitting functions. Since as is now very strongly constrained by measurements 
at LEP and through unpolarised DIS, one might expect that data on gi constrains the polarised gluon densities. 
Unfortunately, errors on gi are large in the low- a; region, where the contribution of the gluons dominate, primarily 
because the asymmetry Ai is small at low-x. For the same reason the flavour singlet sea quark density is also 
rather loosely constrained by data. We investigate the statistics necessary to improve these constraints through DIS 
measurements of at polarised HERA. 

The plan of this paper is the following. In the next section we discuss the various technicalities that distinguish 
different global analyses. This section also serves to set up the notation. This is followed by a section that discusses 
our choice of data used in the fit. The next section contains our results for the LO and NLO fits, and a detailed 
consideration of the parameter errors. A section on some applications of our parametrisations follows this. The final 
section contains a summary of our main results. 

II. CONSTRAINTS ON PARTON DENSITIES 

A. Parton densities and structure functions 

With Nf flavours of quarks, we need to fix 2Nf + 1 parton densities. These are for the 2Nf flavours of quarks 
and anti-quarks and the gluon. For protons or neutrons, the quark and anti-quark densities for the strange and 
heavier flavours are equal. We work with the two (flavour non-singlet) polarised valence quark densities, Vu and Vd, 
corresponding to the up and down flavours. The other non-singlet densities we use are those corresponding to the 
diagonal generators of SU (5) flavour — 

q3 = 2{n-l), qs = 2{n + ~d-2l), 

§15 = 2(S-f J + 35), q2A = 2(Ti + ~d+l + l:~ ■^). (4) 

The initial conditions for evolution arc that below and at each flavour threshold, the density for that flavour of quarks 
is zero. Thus, below the charm threshold we have ^24 — Qi5 — Qo a-nd below the bottom threshold we set §24 = 9o- 
For the singlet quark density, we use 

go^2^|^. =S-^V>, (5) 
/ / 

in preference to the usual S (which is the sum over quark and anti-quark densities of all flavours). The evolution 
equations couple go to the gluon density g. We also define similar unpolarised quark and gluon densities]^. 

Finally, the structure functions a-^d 9ij for the proton and neutron, are given by eq. (H). The unpolarised 
structure functions Ff'" are given by the analogous expression where the polarised parton densities are replaced by 
the unpolarised densities. An isospin flip, interchanging Vu and Vd and switching the sign of gs, relates gf and g". 
After correcting for nuclear effects, the normalised structure function for deuterium is gf = {g^ + gi)/2, and a similar 
expression for Ff. 

We shall have occasion to use the first moments of various polarised parton distributions. We introduce the 
notation — 



^Our convention is that polarised quantities are distinguished from the corresponding unpolarised ones by a tilde. 
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r„(g2)= / dxVuix,Q^), rdiQ^) = [ dxVdix,Q^), 

Jo Jo 

r«(Q')= / dxUx,Q^), rg(Q2)= f dxgix,Q^). (6) 

Jo Jo 

(7) 

We will also use the notation Ty = r„ + and Ty = — Td- The notation 

rr(Q')= C dxg-'^[x,Q^) (8) 
Jo 

is fairly standard. We shall use it in the text. We also use the notation Fa, etc., to denote the first moments of the 
flavoured sea densities. 

B. The fitting strategy 

We use experimental data on the asymmetries A" and Af, measured on proton, neutron ("^He) and deuterium 
targets to constrain the polarised parton densities. We assume full knowledge of the unpolarised parton densities as 
given by some global fit, so that the structure function Fi can be reconstructed using appropriate NLO coefficient 
functions. Then the data on Ai can be converted to gi. We prefer this method to taking the gi values presented 
by experiments, since different experimental groups may make different assumptions about the unpolarised structure 
functions. Such effects would lead to additional normalisation uncertainties in any global fit. 

We have chosen to use the CTEQ4 set of parton densities |^ in our work. We do not expect this choice to affect 
our conclusions strongly since the unpolarised parton densities now have smaller errors than data on the polarisation 
asymmetry Ai. However, with this choice we are constrained to follow some of the assumptions made by the CTEQ 
group — 

1. We work in the MS scheme, since the CTEQ group does that. Other possibilities would have been to work in 
the AB g or JET 0] schemes, but then we would have had to transform the CTEQ distributions. We prefer 
to avoid this procedure, since the best fit parton densities in one scheme do not necessarily transform into the 
best fit densities in another scheme. 

2. We retain the CTEQ choice for the charm quark mass being 1.6 GeV and the bottom quark mass to be 5.0 GeV. 
At each mass threshold, we increase the number of flavours by one, and treat the newly activated flavour as 
massless immediately above the threshold. Parton distributions and as are continuous across these thresholds 

3. We are constrained to use the Aqcd values used in 

4. We take Ql = 2.56 GeV^ in order to avoid having to evolve the unpolarised parton densities downwards. 

In future we plan to study the results of relaxing one or more of these restrictions. 
We follow the parametrisation of CTEQ4 and write — 

fix,Ql)^aax'''{l-xrHl + a3x''^), (9) 

for all densities apart from (73, which is parametrised as 

q3{x,Ql)=aox'''{l-xy'{l + a3V^ + a4x). (10) 

We have made the choice that the large- a: behaviour of any polarised density is the same as that of the unpolarised 
density; in other words, the parameter 02 is the same for corresponding polarised and unpolarised densities (this 
assumption is sometimes given the name "helicity retention property" [|l2|). For simplicity we have also equated the 
polarised and unpolarised values of 04 when this parameter is a power. 

Finally, at Qg we have extended some of the CTEQ assumptions for unpolarised parton densities to polarised. 
These include equating the values of ai for Vu, Vd and 93, taking a4 — I for go, equating the values of a2 for qo and 

qs. We also retain the choice rj = 257 (w + d) — 1/2 in some of our fits, but let it vary in others. Although these 
assumptions seem overly restrictive, the quality of the data does not allow us to fit many of these parameters. We 
discuss some of these points later in this paper. 

The main difference between our parametrisation and previous ones is that we explicitly include a non-zero §3(0;, Qq) 
and break SU{2) flavour symmetry in the polarised sea. This part of the sea density is actually quite well constrained, 
and plays a crucial role in our fits. 
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C. Sum rules 



In a three flavour world, we can write down the following sum rule for the first moments of the nucleon structure 
functions in NLO QCD— 



(11) 



The anomalous dimensions on the right hand side of this equation have been calculated in the MS scheme Q]. Two 
facts used are that in the MS scheme the first moment of the NLO quark coefficient Cg^"* = —2, and the first moment 
of the gluon coefficient function is = 0. The upper and lower signs belong to protons and neutrons, respectively. 
The quantities g^, gs and go are baryonic axial couplings. They are defined as matrix elements of axial vector currents 
between baryon states. Due to the axial anomaly, the singlet axial-vector current is not conserved. As a result, go 
picks up a dependence Hence, go, the moments, and as have to be evaluated at the same in eq. (pT[). 

It is not easy to extract go from low-energy hadron data, although there have been some attempts to do this 
using elastic vp scattering [1^. This gives go — 0.14 ± 0.27. Lattice computations and QCD sum rules p^ ] 
also give similar numbers, but have systematic uncertainties which have to be removed in future. The quantity 
.93 = 1.2670 ± 0.0035 [1^] is obtained from the neutron beta-decay constant. The coupling gs is extracted from the 



decay of strange to non-strange baryons. SU{3) flavour symmetry is used crucially in this extraction |18 . The PDG 
result is qh = 0.579 ± 0.025 ffll. 



Using the coefficient functions in the MS scheme, the moments of the structure functions can be expressed as 



rr(Q') = {i- 



(12) 



where the upper (lower) sign is for the proton (neutron). Apart from eq. (H), we have used the definitions of the 
non-singlet densities in eq. (^ and the singlet parton density in eq. (^. 

There are two sum rules which can be obtained by equating the right-hand sides of eqs. (pi] ) and (|l^). Alternatively, 
we could use some linear combinations. The only one that removes the coupling go is the difference — F", and 
gives the Bjorken sum rule At NLO this is — 

r^ + r3 = g3- (13) 

Note that the first moments of non-singlet densities are independent of order by order to all orders. As a result, 
this equation is also valid in this form to all orders. We impose this form of the Bjorken sum rule on our fits. 

The sum of the two gives the EUis-Jaffe sum rule when the additional assumption g^ = go is made. This cannot be 
correct in QCD because go is dependent and gs is not. Moreover, in the absence of a real measurement of go{Q^), 
an EUis-Jaffe type sum rule cannot constrain the parton densities as the Bjorken sum rule does. Hence, we use such 
a sum rule to extract go rather than to impose it as a constraint on parton densities. 



D. Positivity 

Polarisation asymmetries are the ratios of the difference and sum of physically measurable cross sections. Since 
cross sections are non-negative, asymmetries are bounded by unity in absolute value. In the parton model or in LO 
QCD, these cross sections are directly related to parton densities. Hence positivity of cross sections imply 



/(x,Q2) 



< 1 



(14) 



for the ratio of each polarised and unpolarised density to leading order in QCD. In our LO fits, we impose these 
restrictions. 

At NLO and beyond, this simple relation between parton densities and cross sections no longer holds. Parton 
densities are renormalisation scheme dependent objects; although universal, they are not physical. Hence they need 
not satisfy positivity [l9| , instead one must impose positivity on the actual cross sections. This is numerically difficult 
and requires knowledge of a variety of cross sections evaluated to NLO. Since this knowledge is lacking, and for 
numerical simplicity, we have instead imposed eq. (H) on all our parton density fits. 
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E. Choice of numerical techniques 



Our numerical goal is to evolve parton density functions with absolute errors of at most 10^'^. If this design goal 
were reached, then numerical errors would lie at least an order of magnitude below all other errors. We integrate the 
evolution equations using a 4-th order Runge-Kutta algorithm. The Mellin convolutions required in the evaluation 
of the derivative are computed using a Gauss-Legendre integral. The parton densities are evaluated on a grid and 
interpolated using a cubic spline method. All the numerical algorithms may be found in [po| . 

The knot points of the cubic spline are selected to give an accuracy of 10^^ in the evaluation of the parton densities. 
The Mellin convolutions are also accurate to this order. We require the Runge-Kutta to give us integration errors 
bounded by 10"''. This gives us the error limits we require. We can test these estimates by checking that all sum 
rules are satisfied to within 10"'^. On a 180 MHz RIOOOO processor, the program takes about 0.15 CPU seconds to 
evolve the parton densities by AQ^ = 1 GeV^. 

III. SELECTION OF DATA 

Experiments do not measure the asymmetry Ai directly; they measure the asymmetry between the cross sections 
for lepton and longitudinally polarised hadrons being parallel and anti-parallel — 

_ da(tT)-da(n) 

^ da(n) + rf^(n)' ^ ' 

or a similar asymmetry. At, with transversely polarised hadrons. These asymmetries are related to the two that we 
require by 

Al = D{Ai + r,A2) and At = d{A2 - ^Ai), (16) 

where D and d are depolarisation factors for the virtual photon and ^ and r] are essentially kinematic constants. 
In terms of the ratio of the Compton scattering cross sections for transversely and longitudinally polarised virtual 
photons, 



we can write 



y' + 2[l-y)(l + R) 2-y 



Here 7 — 2Mx/Q <C 1. Using the degree of transverse polarisation of the virtual photon, 
we can write 



(19) 



d = and i = V^- (20) 

Since 7 is very small in the DIS region, the relations Al = DAi and At — dA2 are actually satisfied to high accuracy. 
We then use the further relations 

= (51 - 1^92)/Fi and A2 = 7(31 + 92) /Fi, (21) 

to obtain eq. (Q) when 7^1. It is clear from the second equation that 52 is difficult to measure. 

The main theoretical uncertainty in measurements of Ai is in the values of R used. In fact, many experiments 
use R in two ways. First, it enters the expression for D and cZ, and hence is used to compute QED corrections for 
initial state radiationn and to construct Ai and A2 from Al and At- Next, it is used along with measurements of 



^We thank Abhay Deshpande for drawing our attention to this point. 
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F2 to compute Fi and thus relate gi to Ai. We bypass this second use of R by utihsing experimental data on Ai 
instead of gi. We are forced, however, to accept the first use of i?, In any case, differences between experiments in 
their estimates of D should be factored into the overall normalisation errors. 

The SMC collaboration has data from muon scattering off both proton and deuterium targets. Data was taken in 
separate runs in 1993 and 1996. The most recent publication for Ai is this supersedes previously published data. 
The E-143 experiment at SLAC has data from electron scattering off proton, deuterium and '^He targets. Their most 
recent publication is which supersedes all previous published data on Ai{x, Q^) by this collaboration. The E-154 
experiment at SLAC has data from electron scattering off ^He targets [Esl. The HERMES collaboration in DESY has 



data from positron scattering off protons and He |24 . We have also used data on DIS from ^He taken by the SLAC 



E-142 collaboration [2^]. We have chosen not to utilise data taken by the older EMC collaboration and the E-140 
experiments at SLAC. 

Deuterium is a spin-1 nucleus with the p and n primarily in a relative s-wave state. The d-wave probability is 
estimated to be = 0.05 ± 0.01 p6[ . This is used in the relation between the structure function of deuterium and 
those of p and n — gf = {\ — 3w£)/2)(5^ + 9\)l^- In ^H^, the two protons are essentially paired into a spin singlet, 
and the asymmetry is largely due to the unpaired neutron. Corrections due to other components of the nuclear 
wave- function are small More details are available in pst . 

From the chosen experiments, we have retained only the data on Q^) for > 2.56 GeV^. While this does 

remove some of the low- a; data, the error bars in the removed data are pretty large. We have checked by backward 
evolution that the data which is removed would not have constrained the fits any further. The total number of data 
points used in our analysis is 224. 

In most fitting procedures the statistical errors on measurements are combined in some way with the systematic 
error estimates. Both sets of errors are usually reported in the literature in each bin of data. Whereas this procedure 
is acceptable for statistical errors, it oversimplifies the nature of systematic errors. These latter are correlated from bin 
to bin, and one must use the full covariance matrix of errors in the analysis. In the absence of published information 
on the covariance matrix, one may make the simplifying assumption that the bin-to-bin correlation vanishes, and add 
the statistical and systematic errors in quadrature. This overestimates the errors on data and hence the errors on 
the parameters determined by fitting. We have made a different extremal assumption of neglecting the systematic 
errors altogether. This procedure almost certainly leads us to under-estimate the parameter errors — a point to be 
borne in mind when we discuss large errors and uncertainties in the fits. In summary, our choice of error analysis is 
deliberately conservative. 

Since g2 contains a possible twist-3 contribution, which cannot be written in terms of parton distributions, we 
cannot utilise data on g2 for our fits. However, the twist-2 part is completely determined by gi. In a later section, we 
report an attempt to limit the extent of the twist-3 term using our fitted polarised parton densities. For this we have 
utilised data on proton target from SMC and the E-143 collaboration at SLAC on deuterium target from SMC, 
E-143 and SLAC E-155 Q, and on neutron target from the E-143 and SLAC experiment E-154 [|l|. In all cases, we 
have used the most recent data set and analysis from each collaboration. The quality of data on g2 is poorer than 
that for gi . This is because A2 is small, and extraction of 52 from A2 requires the subtraction of 171 , which itself has 
significant measurement errors. The errors are dominated by statistical uncertainties. 

There remains data from semi-inclusive DIS taken by the SMC ||3^ and HERMES experiments. Analysis of 
these require fragmentation functions and their evolution. Since such analyses are still to reach the stage that 
parton densities have, utilising parametrisations for fragmentation functions would introduce larger errors into our 
fits. For this reason, we have chosen not to use such data in this work. 

IV. RESULTS 

We have made four full analyses — two LO and two NLO, each with and without SU{2) flavour symmetry (denoted 
S and S respectively) for the sea quarks, and with fixed rj = 0.5. In addition, we have made a set of fits with SU{2) 
symmetric sea (§3 = 0) but t] allowed to vary freely. At LO this had no effect on the fit — rj = 0.5 gave the best fit 
and the remaining parameters were identical to LO S. At NLO 77 moved to 0.6, but the parameters remained close to 
the set NLO S. The goodness of fit improved only marginally when 77 was allowed to float. Given the uncertainties in 
the remaining parameters we retain the choice 77 — 0.5 in the main work reported below. 

The goodness of fit, and the constraints imposed by each set of data are summarised in Table R. The values of 
favour the NLO sets slightly. The HERMES proton and the E-154 neutron data (see Figures y, ||) express the 
strongest preferences for the NLO fits; almost the entire change in in going from LO to NLO comes from these 
two data sets. It is probably no coincidence that these two data sets also have the smallest error bars. 

The parameters and their error estimates are shown in Tables It is worth noting that the NLO S densities lie 
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well within the limits allowed by eq. ( p^ , as do all the densities for NLO S except the gluon. The parameters oq and 



as for g in NLO S are at the limit of positivity. The quality of the NLO S fit is shown in Figures |l|-j3| We recommend 
that the NLO parametrisations be used with the CTEQ4M set of unpolarised parton densities and the LO with the 
CTEQ4L set and with appropriate values of Aqcd- 

For the S densities, the normalisation of $3 inherits its error from the valence parameters and the coupling 53. It 
is the best constrained among the parameters describing the sea. Similarly, for the two S densities the normalisation 
of Vd is fixed by the Bjorken sum rule and its error is inherited from the remaining valence parameters. 
In Fig ^ we have shown the variation of 

Ax' =X'- xLn (22) 

when one of the parameters in 5 is varied for fixed values of all the other parameters in the set. The minimum of 
these curves fixes the best-fit value of the parameter, and the points where Ax' = 1 give the 68.3% confidence limits 
on this parameter. Uncertainties in the gluon density are shown in greater detail in Figure ||. This shows contour 
lines of A^^ = 2.3, which encloses the area with 68.3% probability of giving a good description of the data. Also 
shown in the figure are lines of constant Tq^Qq). Although la contours give negative values of F^, we note that 3a 
contours include positive values as well. Note that our error bars are deliberately conservative, due to our neglect of 
systematic errors. Since systematic errors are as large as the statistical errors, the naive procedure of summing them 
in quadrature would have led us to believe that at NLO positive Tg is allowed at 1.5a. 

The huge uncertainties in the gluon distribution due to these parameter variations are illustrated in Figure ^. At 
X = 0.01 the polarised gluon density at NLO can lie anywhere in the range from —10 to —50. This uncertainty at 
low- a; prevents us from investigating this theoretically interesting region. The situation is worse at LO where g is at 
the limit of integrability, since ai ~ —1. Thus F^ is essentially undetermined. 

This large uncertainty in g comes because the only constraint on gluon densities at present are the data on 
variations of gi. Furthermore, the data at a; < 0.1 are most effective in constraining g, and in this range, the data 
have large errors. We have quantified this in Figure ^, where we show selected data on gi in different bins of x as a 
function of (the normalisation has been made arbitrary for ease of viewing) . As the bands of variation due to g 
(in the NLO S set) show, the data does not constrain the gluon density well. 

A few qualitative statements are in order. It is clear that in a scaling theory, DIS data with nucleon or nuclear 
targets can at most fix two linear combinations of quark densities. However, DIS structure functions are Q'^-dependent. 
Hence, in principle, data of arbitrarily high accuracy fixes these two linear combinations at each — i.e., an infinite 
number of functions. In QCD there is a more economical description of this dependence involving the set of 
parton densities given in Section 2. Of course, in the real world data are never infinitely precise, so the question is 
how accurately does evolution fix these densities. Some part of the answer is clear from Figure ^ — improved data 
at low X will constrain g much better and gf, being iso-singlet, would present the best constraint. 

We investigated this question quantitatively by generating fake data at x < 0.1 from the NLO S set. The values 
of so generated were smeared randomly over a 10% band to simulate noise in the data, and error bars of 20-30% 
were assigned to each such data point. This faked set is meant to mimic data that could possibly come from a future 
polarised HERA experiment. We redid our fits with this faked data set replacing all data on for x < 0.1. This 
data brings down the error bars in the parameter ao appearing in 5 by a factor of 4, and improves the errors in oi by 
a factor of 2-3. The error estimates obtained with the fake data set are shown in Figure |[ The region of parameter 
space allowed by this fake data is shown as the grey patch in Figure]^. Taking data at a; < 0.01 or over a larger range 
of Q^, both of which would be feasible at polarised HERA, would constrain g even better [|4j. The ability of RHIC to 
fix the gluon densities is, of course, well appreciated. However, if DIS experiments can fix the gluon densities better, 
then polarised RHIC can be used as a discovery machine. 



V. APPLICATIONS 



A. Flavour asymmetry 

Unlike previous fits of parton densities which had built in the constraint d « m , we have allowed for sea quark 
densities that violate flavour SU (2) symmetry. The fits show that the data tolerate, and even prefer, strong flavour 



symmetry violations. This is easily seen in the first moments of various densities (Table Vl). Since Fq and Fg are 
small, it is clear that the pattern [Fsl <C \Tu\ and Fj « — Fg must follow whenever F3 is large. 

It has been suggested that SU (2) flavour asymmetry in the sea be observed through two combinations of cross 
sections for production of in longitudinally polarised pp scattering — 
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a(w^+,rT)±fT(Ty-,TT) + ^(w^+,n)±^(w^-,n)' ^ ^ 

where A^r (As) is defined with the upper (lower) signs. At LO these asymmetries can be written as the ratio of certain 
combinations of polarised and unpolarised parton densities. At LO the asymmetry Aoy for Drell-Yan pairs can also 
be written in terms of the parton densities. At zero rapidity, Ady is a function of -^r = M/ \/S, where M is the mass 
of the pairs and ^/S is the center of mass energy of the colliding protons. For production at zero rapidity, the 
parton densities have to be evaluated at M = Mw, and hence VS = Mw/Vt. 

In Figure ^ we have shown these asymmetries at zero rapidity as a function of x/S*, or equivalently of y/r. It is clear 
that at \/S appropriate to RHIC, the iso-triplet spin asymmetry, As, is best suited to distinguish the LO S densities 
from LO S. Experimentally studying the dependence of As on over even a limited rang e below y/S ^ 500 GeV 
would be very useful. Ajjy is the least suitable measurement for making this distinction. 

We have already pointed out that the LO fit is unable to decide on the sign of Tg, and that the NLO fits yield 
a negative Tg, although positive values are not ruled out. Although previous fits have seen overlapping ranges of 
allowed Tg, the theoretical bias has been to take large and positive values of Fg. This sign can be easily fixed by 
various experiments at RHIC or in charm production measurements at HERA [ p6[ or the COMPASS experiment in 
CERN I3. 

We would like to caution that parton densities are renormalisation scheme dependent (and hence unphysical) . 
They are universally applicable to all experiments, as long as each experiment is treated in the same scheme [ p8[ . 
Our determination of these densities are in the MS scheme, and statements about their moments are therefore also 
restricted to this scheme. When interpreting the moments of unphysical parton densities, their scheme dependence 
must be held in mind. 



B. Structure functions and couplings 



It is possible to construct physical quantities out of the unphysical first moments of the parton densities. For the 
first moments of the structure functions as given in eq. (p^, wc obtain the values — 

F^" = 0.136 ±0.008, F'/ = -0.053 ± 0.008 (NLO S), 

r{ = 0.157 ± 0.006, F^" = -0.032 ± 0.006 (NLO S) (24) 

at Ql = 2.56 GeV^. These values are within 1.5fT of each other and they compare well with values deduced from 
experiments |j2|,|2|l. 

We can also use eq. (^l]) to extract the value of go. Using as input the above values of F^"'" derived from our fits 
and the PDG value for gs, we find 

,^2n f 0.13 ±0.20 (NLOS), 
^°^^°) = 1 0.35 ±0.15 (NLO si: ^'^^ 

Since go is a physical quantity, it is only to be expected that our determination of go should agree with other analyses, 
such as even if they use some other scheme to arrive at the same result. We will, of course, disagree with them 
on any scheme dependent quantity, such as AT, = F„ + F„ + Fq. Our maximally flavour symmetry violating fits give 
physically reasonable results. 



C. The structure function g2 

Wandzura and Wilczek |3^] have derived a sum rule relating the twist-2 part of g2 to gi — 

52^^(x,Q^) = -51 (-,0^) + ^i^^dy. (26) 

Jx V 

An additional twist-2 contribution to 52, suppressed by the ratio of the quark to the nucleon mass Ml, is ignored 
here. Predictions for the twist-3 contribution have been made using bag models |^l|, QCD sum rules ||42| as well as 
from non-perturbative lattice QCD computations | |4^ . Since some computations predict large twist-3 contributions 
to moments of 52, it becomes interesting to check whether the data on 52 allows such contributions. 
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Figure 10 shows our "prediction" for the twist-2 part of 32 and compares it to measurements of this structure 
function. Clearly the data is compatible with the NLO twist-2 prediction (and also with the parton model result, 
g2 = 0). Between the prediction and the data, there is little room for a twist-3 contribution. Statistics have to be 
improved vastly in order to study higher-twist effects. In fact, COMPASS hopes to make this measurement ^t\ . 

Since the statistical errors are smallest for j it seems that this is the best candidate in which to look for twist-3 
effects. However the data quality needs improvement even here. There is considerable scaling violation in the twist-2 
part of but the large errors prevent any analysis of the Q^-dependence. 



D. The valence densities 



Recently the HERMES collaboration has used semi-inclusive polarised DIS data to extract valence and sea quark 
densities |3^. We have not used these in our fits because this analysis is performed with parton model formalae. 
Nevertheless, it is interesting to compare our fits with these numbers. We display this comparison in Figure |l^. 
The rough agreement is heartening, but the small systematic differences between the fit results and the HERMES 
extraction of the valence densities shows the need for a more accurate QCD analysis of the experimental data, taking 
into account properly the dependence through NLO evolution. 



E. Axion-matter coupling 

We present an example of the application of polarised proton scattering to physics beyond the standard model. 
The Peccei-Quinn solution to the strong CP problem postulates a global symmetry whose spontaneous breakdown 
generates a (nearly) massless pseudo-Goldstone boson called the axion Q]. There is a variant of the original model 
which is still viable [ ^Sf . The axion, a, whose decay constant is fa, couples to fermions, "0/; of mass m/ by the term 

Cint = -igf V'/75V'/ a. (27) 

The coupling gj — Cf{mf / fa). The effective Peccei-Quinn charge, C/ = Xf/N, appears in the coupling instead of 
the actual charge Xj. Here, N is given by '^Xf. There have been several studies Q] of the effective Peccei-Quinn 
charge of the proton. For three flavours, the LO expression can be written as 

/ 

where fJ-f = M/ruf with 1/M = J2^/''^f [0' quarks and leptons is highly model dependent. In the so-called 

KSVZ Q|, and other hadronic axion models, C„ = = Cg = 0. Using quark mass ratios mu/rrid = 0.568 ± 0.042 
and mu/ms — 0.0290 ± 0.0043, obtained by chiral perturbation theory |^^, and our LO S fits, we find that 

Cp = -0.402 (2) and C„ = -0.058 (4). (29) 

The statistical errors in this coupling are dominated by the errors in the fits to polarised parton densities. The 
uncertainty in the quark masses give smaller contributions to these errors. However, the real source of uncertainty 
comes from higher loop corrections. An estimate of this theoretical uncertainty can be obtained by inserting the NLO 
results for various moments into eq. ([2^). This changes Cp by 25% and C„ by 100%. 

The chiral couplings of neutralinos and charginos in generic supersymmetric extensions of the standard model also 
give rise to effective couplings with matter which depend on the moments of the parton densities. Such couplings 
are often needed in astrophysical contexts. Unless these couplings are examined to 2-loop order, they should not be 
evaluated with the NLO moments. 



VI. CONCLUSIONS 



We have made global QCD analyses of data on the asymmetry Ai from polarised DIS without making overly 
restrictive assumptions about the flavour content of the sea. We have extracted polarised parton densities which can 
be used with the CTEQ4 set of unpolarised densities. Our NLO analyses (in the MS scheme) yield the parameter 
sets given in Tables || and HI, and the LO analyses give the sets displayed in Tables IV and 0. Strong SU{2) flavour 
symmetry violation in the sea is supported by the sets of densities in which the first moment of the singlet sea quark 
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density is much smaller than that of the triplet sea quark density. We also found in the NLO fits that the first moment 
of the gluon density is preferentially negative, although positive values are within 2-3 <j of the best fits. The LO fits, 
to the contrary, prefer a positive sign for Vg{Q'^, although negative values are allowed within 2-3 a of the central 
values. All physical quantities obtained from the first moments of our fitted densities have completely sensible values, 
as they must have. 

We have shown that the polarised HERA option may yield highly accurate measurements of the polarised gluon 
densities if measurements of A\ in the range x < 0.125 can be performed with errors of about 25%. Measurements of 
the asymmetry in the iso-triplet part of production at RHIC are likely to be able to pin down the flavour content 
of the sea. Measurements at future facilities for spin physics thus nicely complement each other. 

We have checked that our parametrisations are roughly consistent with semi-inclusive DIS data, although a full 
QCD analysis of this data remains to be performed. We have also shown that these parametrisations, when used to 
determine the twist-2 part of g2 leave very little room for a twist-3 part to this structure function. Finally, we have 
determined the coupling of hadronic axions to matter — an input into several astrophysical constraints on the invisible 
axion. 

In a future publication we plan to make a more detailed study of several issues, including the proper inclusion of 
systematic experimental errors into the analysis and several other technical issues concerning NLO QCD global fits. 

We would like to thank Willy van Neerven for several discussions and clarifications. We also thank the organisers 
of the 6th Workshop on High Energy Physics Phenomenology (WHEPP-6), Chennai, January 2000, where a portion 
of this work was completed. 
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Experiment 


Points 


LO S 


LO S 


NLO S 


NLO S 


SMC (p) 


48 


46.6 


48.4 


43.4 


42.1 


SMC (d) 


53 


53.5 


53.7 


54.3 


50.0 


E-143 (p) 


43 


48.2 


47.8 


48.7 


50.2 


E-143 (d) 


43 


62.1 


60.3 


61.0 


57.0 


HERMES (p) 


9 


23.8 


18.8 


12.6 


15.1 


HERMES (n) 


4 


1.7 


1.6 


1.5 


1.7 


E-142 (n) 


15 


15.9 


16.8 


18.7 


21.1 


E-154 (n) 


8 


22.9 


16.7 


3.3 


2.3 


Total 


224 


274.6 


264.0 


243.7 


239.4 



TABLE L The contribution to ^ from diflerent data sets. The parameter sets marked S impose SU{2) flavour symmetry 
on the sea, whereas the sets marked S do not. 



density 


ao 


ai 


Cl2 


0.3 




Vn 


0.615+1 


-0.32 (±2) 


3.689'' 


12.2 (±2) 


0.873** 


Vd 


-0.61 (±2) 


-0.32" 


4.247'' 


2.2 (±1) 


0.333** 


90 


0.009 (±9) 


-0.2tT 


8.041'' 


8 (±16) 


1.000** 


93 


-0.22^= 


-0.32" 


8.041'' 


0" 


7 (±5) 


g 


-i.otl 


-0.7+1 


4.673'' 


0-2 


1.508* 



TABLE H. The NLO S fits for the parameters in eqs. (1,0) at = 2.56 GeV^ The error estimates shown in the brackets 
apply to the last digit of the estimated value. In case of asymmetric errors, if one of the errors is zero it indicates that the 
parameter is at the limit of positivity. The parameters marked (a) are set equal to some other in the same column, (b) are 
fixed to the value taken by the unpolarised densities, and (c) are fixed by the Bjorken sum rule. 



density 


ao 


ai 


a-i 


as 


Cl4 




1.74 (±1) 


-0.14913 


3.689'' 


3.9ll^ 


0.873** 


Va 


-0.75^= 


-0.149" 


4.247** 


1.6 (±1) 


0.333** 


Qo 


-0.26 (±2) 


-0.08 (±4) 


8.041'' 


6.5 (±7) 


1.000** 


g 


-0.31J 


-o.etl 


4.673'' 


_17+i5 


1.508** 



TABLE HL The NLO S fits for the parameters in eq. (|^ at Qg = 2.56 GeV^. Asymmetric errors and superscripts on the 
numbers have the same meaning as in Table Fll 
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density 




ao 




ai 


"2 




a4 


K 




1.65tg 




-0.159t^ 


3.465'' 




1.146'' 


Vd 








-0.159" 


4.003' 




0.622'' 


Qo 




-0.21 (±2) 




o.oil^ 


6.877'' 




l.OOO'' 


Q3 




o.sr 




-0.159" 


6.877'' 


0" 


-3.5+^ 


9 




-0.16+^ 




-i.otl 


3.666'' 




1.968'' 


TABLE IV. The LO S fits for the parameters in eqs. (|, 0) at C 
numbers have the same meaning as in Table Q 


= 2.56 GeV^ 


Asymmetric errors and superscripts on the 


density 




ao 




ai 


a2 


0.3 








1.91 (±1) 




-0.150 (±3) 


3.465'' 




1.146'' 


Vd 




-1.22^= 




-0.150" 


4.003*" 


0.5tt 


0.622'' 






-0.20 (±2) 




0.03l4 


6.877'' 


2.2« 


l.OOO'' 


9 




-o.mtl 




-i.otl 


3.666'' 


-I5tf 


1.968'' 


TABLE V. The LO S fits for the parameters 
numbers have the same meaning as in Table Q 


in eq. (||) at Qq 


= 2.56 GeV^ 


Asymmetric errors and superscripts on the 






LO S 




LO S 




NLO S 


NLO S 


r„ 




0.875 (±5) 




0.829 (±1) 




0.909 (±8) 


0.85 (±3) 


Td 




-0.40 (±4) 




-0.338 (±6) 




-0.36 (±1) 


-0.52 (±3) 


To 




-0.029 (±3) 




-0.028 (±3) 




-0.058 (±6) 


-0.003 (±3) 














-0.2tl 


-1.6 ± 1.0 


La 








0.107 (±3) 






-0.10 (±2) 






-0.0059 (±6) 




0.0212 (±9) 




-0.012 (±1) 


-0.024 (±6) 






-0.0059 (±6) 




-0.0323 (±9) 




-0.012 (±1) 


0.024 (±6) 


2r,- 




-0.0029 (±3) 




-0.0027 (±3) 




-0.0058 (±6) 


-0.0003 (±3) 



TABLE VI. Moments of various densities at Qq — 2.56 GeV^. By our initial conditions Fg — 2ro/5. The numbers in 
brackets are the errors on the last digit of the central value. Fg is essentially undetermined at LO. 
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FIG. 1. Data and fits for A^". The data are at different Q^, but the curve is the asymmetry from NLO S set calculated at 
fixed Ql = 5 GeV^ 



+ HERMES 
□ E-142 
.E-154 










* ' 1 — 











T5l3T 0!l 



X 

FIG. 2. Data and fits for A". The data are at different Q^, but the curve is the asymmetry from NLO S set calculated at 
fixed Ql = 5 GeV^ 




FIG. 3. Data and fits for Af. The data are at different Q^, but the curve is the asymmetry from NLO S set calculated at 
fixed Qg = 5 GeV^ 
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FIG. 4. Plots of Ax^ against the parameters ao and ai for g in the two NLO sets. In the set NLO S, ao is at the boundary 
of positivity. 
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FIG. 5. The covariance of the fitted parameters ao and ai for g for the two NLO sets with as kept at their respective best fit 
values. The crosses show the best fit points, and full lines are the contours enclosing the 68% confidence limits. The reference 
points A and B are used to quantify the variation in gluon densities in Figure |^. Along the dashed lines Fg — —1/2", for the 
values of n marked. The grey patch is the region allowed by the faked data discussed later. 
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FIG. 6. The absolute values of the polarised valence quark densities are shown along with the range of allowed \g\ in the 
NLO S set. The gray band showing this uncertainty is the band enclosed by the densities obtained at points A and B marked 
in Figure H. 
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FIG. 7. We show gi{x,Q^) in selected bins of a; as a function of log(Q^/A^). The data on from SMC are shown by open 
circles and from E-143 by filled circles; on gf from the same experiments by open and filled triangles respectively, and on gf 
from E-142 by open squares and from E-154 by filled squares. The bands are the uncertainty in gi induced by the uncertainty 
in g shown in Figure m. Data in different bins of x are offset vertically for visibility. 
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FIG. 8. Plot of Ax against the deviation from the best fit value of the gluon parameters ao and ai using faked data with 
20-30% errors in measurements of Ai for x < 0.1. 
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FIG. 9. The asymmetries Aby and Aa and Af, (defined in eq. ^3|) computed at LO. The full lines are obtained with LO S 
and the dotted lines with LO S. The bottom scale is for A^^^ as a function of \fS and that at the top for Aby as a functions 
of The vertical band marks out the range 200 GeV < \/S < 500 GeV. 
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FIG. 10. Data on the structure function compared with the twist-2 predictions of eq. I 
parametrisation evolved to = 5 GeV^. From left to right, the figures are for and g^. 



evaluated with our NLO S 




or 



FIG. 11. The asymmetries (u + it)/(ii + it) (boxes) and {d -\- d) j id -\- d) (pluses) extracted by a parton model analysis of 
experimental data |33] (the two overlapping points at a; = 0.35 have been separated for clarity) compared to our NLO fits at 
Q2 = 5 GeV^. NLO S is the full line and NLO S is the dotted line. 
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